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The S-matrix theory formulation of closed-orbit theory recently proposed by Granger and Greene 
is extended to atoms in crossed electric and magnetic fields. We then present a semiclassical quan- 
tization of the hydrogen atom in crossed fields, which succeeds in resolving individual lines in the 
spectrum, but is restricted to the strongest lines of each n-manifold. By means of a detailed semi- 
classical analysis of the quantum spectrum, we demonstrate that it is the abundance of bifurcations 
of closed orbits that precludes the resolution of finer details. They necessitate the inclusion of uni- 
form semiclassical approximations into the quantization process. Uniform approximations for the 
generic types of closed-orbit bifurcation are derived, and a general method for including them in a 
high-resolution semiclassical quantization is devised. 
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I. INTRODUCTION 

Closed-orbit theory was first introduced by Du and 
Delos § and Bogomolny Q some twenty years ago to 
interpret the modulations observed in the photo-absorp- 
tion spectra of hydrogenic Rydberg atoms in a magnetic 
field close to the ionization threshold. Since that time, 
it turned out to be a powerful and flexible tool for the 
semiclassical interpretation of a variety of spectra. It 
has been used to describe atoms in electric j| as well as 
parallel j|, ^ or crossed M, [?], |sj electric and magnetic 
fields. In the case of non-hydrogenic atoms, the influ- 
ence of the ionic core can be modelled either by means 
of an effective classical potential |], |l0| or in terms of 
quantum defects |ll], |l 2fj . Recently, closed-orbit theory 
has even been shown to be applicable to the spectra of 
simple molecules in external fields Jr| . 

A complete description of photo-absorption spectra re- 
quires the calculation of the energies E n of the excited 
atomic states and the strengths of the spectral lines, 
which is characterized by the dipole matrix elements 
(i\D\n) between the initial state \i) and the Rydberg 
state |n), where D is the component of the dipole opera- 
tor describing the polarization of the exciting laser field. 
These quantities are neatly summarized in the response 
function 
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denotes the retarded Green's function. 

Closed-orbit theory provides a semiclassical approx- 
imation to the quantum response function ([!]), which 
splits into a smooth part and an oscillatory part of the 
form 
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where the sum extends over all classical closed orbits 
starting from the nucleus and returning to it after having 
been deflected by the external fields, S c . . is the classi- 
cal action of the closed orbit, and the amplitude A c .o. 
describes its stability and its starting and returning di- 
rections. Its precise form depends on the geometry of 
the external fields. In section III, it will be specified for 



systems with and without a rotational symmetry. 

Although the closed-orbit sum (||) appears to provide 
a straightforward means of calculating the response func- 
tion from the classical closed orbits, this is actually not 
the case because the sum usually diverges for real en- 
ergies E. Thus, the quantal information cannot be ex- 
tracted directly from the semiclassical expansion. One 
particular and widely applicable method to overcome the 
convergence problems of the closed-orbit sum is semiclas- 
sical quantization by harmonic inversion Jl4], ^5). For 
the hydrogen atom in a magnetic field, this method has 
been shown |l(| to be capable of extracting semiclassi- 
cal eigenenergies and transition matrix elements from a 
closed-orbit sum. 

In the present paper we will investigate how these re- 
sults can be generalized to the hydrogen atom in crossed 
electric and magnetic fields. This problem is considerably 
harder than the treatment of the diamagnetic hydrogen 
atom, which possesses a rotational symmetry around the 
field axis. Due to that symmetry, in classical mechan- 
ics the angular momentum around the field axis is con- 
served. So is, in quantum mechanics, the magnetic quan- 
tum number m. In crossed fields, the rotational symme- 
try is broken. As a consequence, the selection rules for 
the m-quantum number no longer hold, and a multitude 
of additional lines appears in the quantum spectrum. At 
the same time, the determination of classical closed or- 
bits gets significantly more difficult because three non- 
separable degrees of freedom have to be dealt with. A 
detailed description of the intricate pattern of closed or- 
bits and their bifurcations was given in an accompanying 
paper |T7f| . That data forms the basis of the present 
work, where the semiclassical treatment of the crossed- 
fields hydrogen atom will be dealt with, and we will freely 
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use the nomenclature introduced in [ fL7| . 

After the essential properties of the crossed-fields 
Hamiltonian have been summarized in section 0, wc 
start, in section III, with a derivation of the closed-orbit 
formula (||) in the context of the S'-matrix formulation of 
closed-orbit theory introduced recently by Granger and 
Greene We show that the novel framework can 

be extended to the crossed-fields situation, and we clar- 
ify some misleading conclusions arrived at in jl8). Sec- 
tion IV describes the quantum spectrum under study, 
and section ^ compares it to a semiclassical spectrum 
in low resolution. In section VI, the results of a high- 
resolution semiclassical quantization using the technique 
of harmonic inversion are presented. The semiclassical 
spectrum correctly identifies the strongest spectral lines, 
but it fails to describe finer details of the quantum spec- 
trum. In section VII, wc compare a quantum recurrence 



spectrum to the classical data to show that the prin- 
cipal source of this difficulty lies in the abundance of 
closed-orbit bifurcations. Uniform approximations pro- 
vide a tool to cope with the divergences introduced into 
semiclassical spectra by bifurcations of classical orbits. 
A general technique for their construction is described 
in section VIII, and uniform approximations for the two 



types of generic codimension-one bifurcations identified 
in |17t are derived. Finally, in section IX we demonstrate 



how uniform approximations can be incorporated into 
recurrence spectra, thus paving the way for their inclu- 
sion into the high-resolution semiclassical quantization 
by harmonic inversion. 



II. THE CLASSICAL HAMILTONIAN 

Throughout this work, we will assume the magnetic 
field to be directed along the z-axis and the electric field 
to be directed along the x-axis. In atomic units, the 
Hamiltonian describing the motion of the atomic electron 
then reads 
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where B and F denote the magnetic and electric field 
strengths, respectively, r 2 = x 2 +y 2 +z 2 , p 2 = x 2 +y 2 , and 
L z is the z-component of the angular momentum vector. 
By virtue of the scaling properties of the Hamiltonian (|4|) , 
if all classical quantities are multiplied by suitable powers 
of the scaling parameter 



(5) 



the dynamics can be shown not to depend on the energy 
E and the field strengths B and F separately, but only 
on the scaled energy E = w 2 E and the scaled electric 
field strength F = w 4 F. In particular, classical actions 
scale according to S = wS. Thus, the semiclassical limit 
of large classical actions corresponds to the limit of large 
w. 



The way of recording a quantum spectrum which 
is best suited for semiclassical investigations is scaled- 
energy spectroscopy. A spectrum then consists of a list 
of the scaling parameters w n characterizing the quantum 
states for given scaled energy E and scaled electric field 
strength F. Scaled-energy spectroscopy offers the ad- 
vantage that the underlying classical dynamics does not 
change across the spectrum. It will be adopted through- 
out this work. 



III. THE S'-MATRIX FORMULATION OF 
CLOSED-ORBIT THEORY 

A. General formalism 

The basic observation fundamental to all of closed- 
orbit theory is a partition of space into physically distinct 
regions. In the core region close to the nucleus, the Ryd- 
berg electron interacts in a complicated manner with all 
electrons of the ionic core. This interaction is manifestly 
quantum mechanical in nature, it cannot be described 
in the framework of semiclassical theories. On the other 
hand, the interaction of the Rydberg electron with the 
external fields is much weaker in the core region than its 
interaction with the core, so that the fields can safely be 
neglected. Therefore, a description of the core obtained 
in the field-free case can be used. In particular, the ini- 
tial state of the photo-absorption process is assumed to 
be localized in the core region and not to be influenced 
by the external fields. 

In the long-range region far away from the nucleus, on 
the other hand, the external fields play a dominant role, 
whereas there is no interaction with the ionic core ex- 
cept for the Coulomb attraction of its residual charge. In 
this region, the dynamics of the Rydberg electron is well- 
suited for a semiclassical description. It is independent 
of the details of the ionic core. 

In order to establish a link between the dynamics in 
the core and long-range regions, a matching region is as- 
sumed to exist at intermediate distances from the nucleus 
where both the external fields and the interaction with 
the core are negligible. Thus, in the matching region 
the simple physics of an electron subject to the residual 
Coulomb field of the core is observed. 

Recently, Granger and Greene [|l8| proposed a novel 
formulation of the theory based on ideas borrowed from 
quantum-defect theory. Their formulation achieves a 
clear separation between properties of the external field 
configuration and the ionic core, which are encoded in 
separate S'-matrices. Suitable approximations to the core 
and the long-range S'-matrices can be derived indepen- 
dently. Therefore, the formalism can be expected to al- 
low a generalization of closed-orbit theory to atoms with 
ionic cores exhibiting more complicated internal dynam- 
ics than have been treated so far. 

The derivation given by Granger and Greene treated 
the case of an atom in a magnetic field only. It will now 
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be extended in such a way that it holds for combined 
electric and magnetic fields with arbitrary field configu- 
rations. To this end, the ansatz and basic formulae of 
Granger and Greene's theory will be summarized in this 
section. A more detailed treatment can be found in their 
paper jjq| . In subsequent sections, we will then turn to a 
discussion of the long-range scattering matrices pertinent 
to different external field configurations. 

To lay the foundation for a definition of the S-matrices, 
we pick a basis set v[<™ rc and \1/^ R of wave functions of 
the Rydberg electron valid in the core and long-range 
regions, respectively, and expand in terms of spherical 
harmonics 

* c fc ore (LR) (r, *,<p) = \Y, Yk ' (*> ^K'T (LR) (r) ■ (6) 

k' 

The channel index k is to be understood as a double in- 
dex (I, m) characterizing the spherical harmonics. When 
studying a complicated atom with more than one relevant 
state of the core, additional information can be included 
in the channel functions Y k . 

In the matching region, the radial function matrices 
F COIC and _F LR can both be expressed in terms of ra- 
dial Coulomb functions. We use the functions f k {r) and 
fk( r ) satisfying outgoing and incoming wave boundary 
conditions, respectively, given b y ||19| and choose the ra- 
dial functions to be of the form ||4.5| 

F™*(r) = i [f+(r) S c k ?r - / fe "(r) 5 Wk ] , (7) 

F^(r) = -i [#(r) 5 k , k - /£(r) ■ (8) 

Physically, these choices mean that the basis function 
,j,core j g a superposition of a single incoming wave in chan- 
nel k and the outgoing waves in different channels pro- 
duced from it by scattering off the core. Similarly, \&j£ R 
consists of an outgoing wave in channel k and the return- 
ing waves generated by scattering off the external fields. 
The scattering matrices S_ corc and 5 LR thus summarize 
the physical properties of the core and the external fields, 
respectively. They are determined by the condition that 
the radial functions obey suitable boundary conditions, 
i.e. F_ coic is regular at the origin, whereas vanishes 
or satisfies outgoing-wave boundary conditions at infin- 
ity for bound and free states, respectively. For hydrogen, 
S_ colc is the identity matrix. 

Following previous work by Robicheaux |20| , Granger 
and Greene derive the following expression for the re- 
sponse function ([!]): 

9 = id} ( 1 + 2 (5 COTC 5 LR ) + 2 (^S cole S LR ) 2 

+ 2 (S coie S hK ) +...)d, 

where the vector d comprises the energy-dependent 
dipole matrix elements 

d k (E) = (*~»(25)|Z>|t) (10) 



between the initial state and the core-region channel wave 
functions. For hydrogen they can be computed explicitly 
(see, e.g., § or [[§). 

The terms of the series (|9|) embody contributions from 
paths where the Rydberg electron takes zero, one, two, 
etc. trips out into the long-range region and back to the 
core before interfering with the initial outgoing wave. In 
the semiclassical approximation, £ LR will be given in 
terms of closed orbits. A returning wave is associated 
with each returning classical orbit. By a general ionic 
core, it is scattered into all directions. The parts of the 
wave scattered into the outgoing direction of a closed or- 
bit will then follow this orbit until they return to the 
core again. Thus, core scattering leads to a concatena- 
tion of different closed orbits Jll], [l2|. In hydrogen, the 
Coulomb center scatters the incoming wave back into its 
direction of incidence, so that there is no coupling of 
closed orbits. Terms describing repeated scattering off 
the external fields are therefore absent from the sum, 
and the hydrogen response function can be decomposed 
into a smooth part 

g = id?d, (11) 

which is the same as in the field-free case and contains 
"direct" contributions where the electron does not scatter 
off the external fields at all, and an oscillatory part 

<? osc = 2i J S LR d (12) 

generated by the electron going out into the long-range 
region and being scattered back to the nucleus. It is this 
part which describes the impact of the external fields. 

The basis for a semiclassical approximation is provided 
by the retarded Green's function G(x, x'\ E) describing 
the propagation of the electron from x' to x at the energy 
E. It can be expanded in terms of the channel functions 
as 

G(x, x'; E) = -L Y k (0, <p) G kk , (r, r';E) Y£(#, J) 

kk' 

(13) 

with 

G kk , (r, r'; E) - rr' (k\G(x, x'; E)\k') . (14) 

The long-range scattering matrix is related to the Green's 
function matrix by |l8fl 

s LR = -[rMr'G^roMrMr 1 , (15) 

tir — — 

where ro is the matching radius, /~ is the diagonal ma- 
trix 

f kk >(r)=f k (r)Skk> (16) 

comprising the radial wave functions, and GJr, r') de- 
notes the part of G(r, r') satisfying incoming-wave 
boundary conditions at the final radius r. The latter con- 
dition ensures that only electron paths approaching the 
matching radius from the long-range region contribute 
to 5 LR , whereas paths that traverse the core region are 
omitted. 
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B. Closed-orbit theory for crossed-fields systems 

To obtain a semiclassical approximation to the long- 
range scattering matrix, we make use of the semiclassical 



Green's function derived by Gutzwiller 22 
2tt 



G sc \x 7 x';E) = 




D\ exp (iS 
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l 2 a 



(17) 



where the sum extends over all classical trajectories lead- 
ing from x' to x at the energy E, n is the number of 
degrees of freedom, S is the classical action along the 
trajectory, a the number of caustics along the trajectory, 
and 
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d 2 S 


dxdx' 
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(18) 



is the amplitude for the contribution of the trajectory. 
By (|l4|), we obtain a semiclassical approximation to the 
Green's function matrix 



Gf£,(r ,r ;E) 



2tt 



(2ttz) 2 ° 

xy fe *(^)*V0?V) 



class, traj 
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(19) 



As usual in semiclassics, the integrals will be evaluated 
in the stationary-phase approximation. It yields a sum 
over all classical trajectories leaving the matching sphere 
at a direction given by (fiijipi) and returning to it at 
(•df^ipf). The condition that G(ro,ro) obeys incoming- 
wave boundary conditions at the final radius translates 
into the condition that only orbits going out from the 
matching sphere into the long-range region and then re- 
turning to ro are to be included, whereas orbits passing 
through the core region are omitted. If all factors in the 
integrand except for the exponential are assumed to vary 
slowly, the stationary-phase approximation reads 

G s k t (r Ql r ;E) = 2irr 2 ^ sin ^ sin f 



det 



d 2 S 



x exp 



/ 7T 

\iS(r ,r ) - i-(a + k) 



(20) 



where k is the number of negative eigenvalues of the Hes- 
sian matrix of S occurring in the prefactor. 

Because the initial state is assumed to be well local- 
ized, it is clear that the outgoing waves generated by 



the photo-excitation originate in the immediate neighbor- 
hood of the nucleus. Therefore, only trajectories leaving 
the matching sphere radially need to be included in fl2C| ) . 
By the same token, the trajectories can be assumed to 
return to the matching radius radially. Thus, they are 
parts of closed orbits starting precisely at the nucleus 
and returning there. 

By transforming jl8| ) to spherical coordinates and 
making use of the relations 



dS_ 
dx 



V 



dS_ 
~dE 



f 



(21) 



the amplitude factor D for radial trajectories can be sim- 
plified to 



D = — 



1 



det 



rr' r 2 r' 2 sin § sin §' d("d,(p) 
The determinants occurring in (|2(i| ) combine to 

d(p'#,p') f, . d 2 S 



(22) 



det ■ 
= det 
= det 



dct 



d(p'#,p'<p,P0,P<p) 
d(fi,<p,p-d,P<p) 

d(p#,p v ) ' 



det 



9(-p'^,-p' v ,P-a,P v ) 



(23) 



With these results, the Green's function matrix assumes 
the form 



■y/sin Wi sini?/ 



iS {vq ,r$) —i-n {a -\- k) / 2 

(24) 

X = fc- 



det 



d(po,,P Vf ) 



The determinant in the denominator of (|24j) measures 
the dependence of the final angular momenta of the tra- 
jectory upon the starting angles. As it stands, it suffers 
from the singularities of the spherical coordinate chart: 
At the poles, neither the angle (p nor the angular mo- 
menta p§ and p v are well defined, so that close to the 
poles, the calculation of the determinant becomes nu- 
merically unstable. The determinant can be rewritten in 



the form 23 



dct 



d(Pj} f ,P Vf ) 



= sin "di sin flf M 



(25) 



with a 2 x 2-determinant M devoid of any singularities. 
The parameter M was already used in jlTj to study the 
bifurcations of closed orbits. We showed there that a 
closed orbit bifurcates if and only if M = 0. With the 



■5 



form (|25|) of the stability determinant, the semiclassical 
Green's function matrix reads 



G 



x exp 



(iS{ra,ra) - i^(a + k)) , 



(26) 



which is free of any singularities introduced by the spher- 
ical coordinates. 

By virtue of (|l^) , the semiclassical long-range scatter- 
ing matrix reads 



Skk' — %i > ■ 



1 



1 



1 



Yj;('& f ,(p f )Y kl ('& i ,(p i ) iS(ro , ro) _ i7r(CT+K) / 2 

A , fc- 

^/\M\ 



(27) 



This expression can be further simplified if, for excited 
states close to the ionization threshold, the radial wave 
functions 



/f(r) « -iVtH$UV8r) 



(28) 



are approximated by the zero-energy wave functions, and 
the Hankel functions are replaced with their asymptotic 
forms for large arguments p4| 



2 / 7T 

exp I — ix + i—v 



(29) 



This approximation has proven accurate in many cases of 
interest, but it was called into question by Grange r and 
Greene [Q~8| . It will be discussed further in section III C , 



where we will show that there is no reason to doubt its 
reliability. It leads to 



cLR 



CO. 



1) i+i' YZ n (0f,<Pf)Yi> m >(#i,<Pi) 



x exp 



(i (S(r , r ) + 2V8n7) - i|(<r + «)) 



(30) 



because, due to the conservation of energy, f 2 /2 = 1/r if 
E = 0. In equation ^0|, the channel indices k = (I, m) are 
finally written out explicitly. 

For a radial trajectory in a hydrogen atom going out 
from the nucleus to r = tq at zero energy, the action is 
V8r , so that 



S c .o. = S(r , r ) + 2y/8r^ 



(31) 



is the action of a closed orbit, measured from its start at 
the nucleus to its return. The semiclassical long-range 
S-matrix finally reads 



LR 



-2ttV(-1) 



CO. 



x exp 



(iS c .o. - v^ip + 



(32) 



Both the action Sc.o. and the stability determinant M 
are evaluated at the nucleus rather than on the matching 
sphere. The response function is given by 



(33) 



x exp ( iSc.o. - i—fJ> 



where the Maslov index fi = a + k + 1 was increased by 
1 to absorb an additional phase, and the function 



(34) 



with the core-region matrix elements di m given by (|1C|), 
characterizes the initial state and the exciting pho- 
ton. Through the di m , the function y(d,<p) is energy- 
dependent. In accordance with the choice of zero-energy 
radial wave functions in the S'-matrix elements, y("d, if) 
will be evaluated at zero energy. This approximation has 
proven accurate in all applications of closed-orbit theory 
considered in the literature so far. However, from the 
S'-matrix theory derivation it is obvious that the energy- 
dependence of both the dipole matrix elements di m and 
the S-matrix elements can easily be included should the 
need arise. The semiclassical response function (|3^) has 
the anticipated form (|^) with 



•Ac 



An 



(35) 



C. Closed-orbit theory for symmetric systems 

An atom in a single (electric or magnetic) external field 
possesses a rotational symmetry around the field axis, 
which must be taken into account in the derivation of 
the closed-orbit formulae. The symmetry gives rise to 
a conserved magnetic quantum number to, so that the 
angular momentum quantum number I remains the only 
relevant channel index. The semiclassical scattering ma- 
trix reads Oil 



S^ R = 2 3 /V/ 2 



y/\A\ sin i3i sin?9/ 



X ^(^,0)^,(^,0) e 'S('-o ! ro)- lI r,V2-3 M /4 



where 



(36) 



(37) 



jl is the number of poles of A encountered along the 
trajectory, and the sum includes all classical trajectories 
with azimuthal angular momentum m joining the circles 
given by polar angles $j and df on the matching sphere. 
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If the radius of the matching sphere is much larger than 
the extent of the initial state, the trajectories can again 
be assumed to leave the sphere and return to it radi- 
ally. Strictly speaking, this condition can only be met if 
m = 0, which we will assume in what follows. If m ^ 0, 
the initial angular velocity ip must be non-zero, but it 
will be small if the matching radius is large. In this case, 
the trajectory will not actually close at the nucleus, but 
swing by at a short distance. 

Using, as above, the radial wave functions at zero en- 
ergy, we obtain the semiclassical scattering matrix 



^ R = -(2 7 r) 3 / 2 (-l) w ' l ^ v /|A|sin^sin,? / 

™. V (38) 

X Y? m {df,Q)Y Vm >(0i,O) e «c.o.-ftr/i/2-3iir/4 

and the response function 



g osc (E) = 2(2nf 2 Y^\/\ A \ sin ^ sin ^/ W/Wi) 

CO. 

x exp (iS c . . - + i^j 



with fj, = p, + 2 and 

W = £(-l)'di*W0,O) 
i 

This result has the form (|J) with 



(39) 



(40) 



A.o. = 2(2tt) 3 / 2 ^/| .4| sin ^ sin / ^* (0 f )y (tf* 

0? 



x exp 



(41) 



It differs from the result obtained previously by Du and 
Delos [Q in that in their work the amplitude factor A of 
( |37| ) is replaced with 



.4 




(42) 



This discrepancy was noted and numerically investigated 
by Granger and Greene . They attribute it to the ap- 
proximation of using zero-energy wave functions, which 
can easily be avoided in the iS-matrix theory, but is an 
integral part of the derivation given by Du and Delos. 

For the closed orbit perpendicular to the field in the 
diamagnetic Kepler problem and a scaled matching ra- 
dius of fo = 0.01, the amplitudes d37j ) and ( ^ ) are plot- 
ted in figure [[} This figure is similar to figure 1 in [5j| , 
although for the latter the matching radius is not given. 
The agreement is excellent at scaled energies close to 
zero, but becomes poor if the energy decreases. How- 
ever, contrary to their conclusions, the lack of agreement 
is not due to the zero-energy approximation, but rather 
to the dependence of the amplitudes on the matching 
radius. 
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FIG. 1: Scaled semiclassical amplitude factors after Granger 
and Greene |L8| (Equation (|37|), solid line) and after Du and 
Delos Jl] (Equation (^), dashed line) for the closed orbit 
perpendicular to the magnetic field as a function of the scaled 
energy. Fhe scaled matching radius is fo = 0.01. 



This statement can be verified most conveniently if the 
motion is described in semiparabolic coordinates 



[i = \fr + z , v = \pr- 



(43) 



If the trajectory is recorded as a function of a parameter 
r related to the time t by 



dt = 2r dr 



(44) 



and a prime denotes differentiation with respect to r, for 
trajectories with vanishing azimuthal angular momentum 
the equations of motion in the Coulomb region read 



p' = 2Efi , p' v = 2Ev 



(45) 



These equations are devoid of any singularities, so that 
they can conveniently be used to discuss the motion close 
to the nucleus. The transformation inverse to ([l3|) is 
given by 



u 2 - v 2 

■d — arccos — . (46) 



/i 2 + v 2 



The momenta transform according to 



Pr 



fiPft + Vp v 
/i 2 + V 2 



P-a 



[ip v - vp^ 
2 sign(/iz/) 



(47) 



Note that the transformation from semiparabolic to 
Cartesian coordinates is not one-to-one, but that \x and 
v are fixed up to the choice of sign only. 

To evaluate ( |37| ) and (^|), the derivatives dp^Jd-di 
and dflf/ddi, must be calculated and their dependence 
on the matching radius r must be determined. As the 
radial trajectory specified by a starting angle is inde- 
pendent of the radius where the angle is measured, the 
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r-dependence of the derivatives is determined by the re- 
turning trajectories only. It can be evaluated as follows: 
We arbitrarily fix the returning time of a closed orbit 
at r = 0, so that /x(0) = v(0) = 0. The solution to (gf) 
describing a trajectory returning at an angle is given 
by 



2 



M(r) 


_ n cos(i? / /2) 


sin 


u(t) 


_ n sin(^/2) 
V-2E 


sin 


pA t ) 


&f 

= 2 cos cos 






„ . ^ 
= 2 sin — cos 
2 





-2Et 



r— . df 

'2rsin — 
2 



-2£r 



2-y/T 



2VT 



- Er cos ^ 

0f 

£r cos — 
2 



(48) 

where the coefficients were chosen to satisfy the conser- 
vation of energy and to give the correct returning angle 
after a transformation to Cartesian coordinates. The sec- 
ond expression in each line follows from /x 2 + v 
whence for r < 



2r, 



sin (V-2Et\ = -\/-Er ., 
cos (v^rj = v 7 ! + Er 



(49) 



Equations of motion for the derivatives d/i/dfli and 
dv/dfli are obtained by linearizing (f45|). Since Q4q ) is al- 
ready linear, the derivatives satisfy the same equations of 
motion as the coordinates themselves as long as the elec- 
tron moves in the Coulomb region. There the solutions 
read 



-f = ^4=sin [V-2Et +^cos [V^ZEt) (50) 



and 



= a M cos (y-2E t) - V-2E 6 M sin (V-2E r 



(51) 



Equation (|g) yields 



+ b ll VTTEr , 



% = a (1 x/iT& - V^Eb, 



(52) 



so that the coefficients 



a u = 



(53) 



can be identified with the values of the derivatives ob- 
tained at r = 0. Analogous expressions hold for dv/dfii. 
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dp-9 



1 



2sign(/i^) \^ d'di 

1 (m„ 



2 sign(/w) \ <9#i 



dv 



(54) 



can be evaluated. It is independent of r, as could have 
been anticipated from the fact that p$ is a component of 
the total angular momentum and thus is conserved along 
the trajectory once the electron has entered the Coulomb 
region. The amplitude A~ l can also, up to an immaterial 
choice of sign, be identified with the monodromy matrix 
element 



1 / dvt d[it 
mi2= 2\W^-m V »> 



(55) 



introduced by Bogomolny || to describe the semiclassical 
amplitudes, so that the amplitudes derived by Granger 
and Greene from the 5*-matrix theory agree with Bogo- 
molny's. 

Similarly, the amplitude (^) used by Du and Delos 
reads, by ©, 



1 



r dd 
2~Wi 
sign(^) 



r dpu f 



-Pu f 



d0 
r dp llf 



- Vl + Er 



Off 



- y/l + Er 



(56) 



Thus, the amplitudes A and A\ agree in the limit of 
vanishing matching radius, but the amplitude A\ pro- 
posed by Du and Delos exhibits a strong dependence on 
r, whereas the amplitude A given by Granger and Greene 
does not. These findings can also be confirmed numeri- 
cally. Figure || shows the two amplitudes for the closed 
orbit perpendicular to the magnetic field at a scaled en- 
ergy of E = —2 as a function of the scaled matching 
radius fo. The dependence of A\ on fo is considerable. 

We have thus shown that, contrary to the conclusion 
reached by Granger and Greene, the discrepancy be- 
tween their semiclassical amplitude and that obtained 
by Du and Delos is not due to the zero-energy approxi- 
mation, but rather due to the choice of a finite matching 
radius. In addition, the amplitude derived by Granger 
and Greene is not specific to the S'-matrix formulation, 
it agrees with the result derived earlier by Bogomolny in 
the context of a semiclassical wave function formalism. 
Nevertheless, as it eliminates the need to specify a finite 
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FIG. 2: Scaled semiclassical amplitude factors after Granger 
and Greene (Equation (|37[), solid line) and after Du and 
Delos jjj (Equatiion (^) , dashed line) for the closed orbit per- 
pendicular to the magnetic field as a function of the matching 
radius at E = — 2. 



matching radius and allows one to calculate all classical 
quantities at the nucleus, it seems more appropriate than 
the amplitude given by Du and Delos, which introduces 
a certain arbitrariness in the choice of a matching radius. 



IV. THE SCALED QUANTUM SPECTRUM 

If Schrodinger's equation for the crossed-fields hydro- 
gen atom is rewritten in terms of the scaled energy and 
the scaled electric field strength, a quadratic eigenvalue 
problem for the scaling parameter w is obtained. An ex- 
act numerical method of solution for the quadratic eigen- 
value problem has become available only recently p5[ . 
We resort to the method introduced by Main Q] , which 
relies on an approximate linearization of the eigenvalue 
problem to compute eigenvalues in a small spectral in- 
terval. The accuracy of the linearization can be verified 
by comparing results calculated using different overlap- 
ping intervals. The eigenvalues are obtained to a relative 
accuracy of at least 1CP 7 , which is far beyond the typi- 
cal accuracy of semiclassical approximations, so that the 
algorithm is well suited to this study. 

In the following we will discuss quantum and semi- 
classical photo-absorption spectra obtained for the scaled 
energy E — —1.4 and the scaled electric field strength 
F = 0.1 with the initial state |2p0) and light linearly po- 
larized along the magnetic field axis. A quantum spec- 
trum for these parameter values is show n in figure ^[ As 
for a semiclassical analysis (see section VII) it is essen- 
tial to have as many eigenvalues available as possible, the 
calculation was extended up to w — 100. The spectrum 
shown in figure [| contains nearly 30,000 lines, many of 
which are too weak to be discernible in the plot. 

The eigenenergies of the field-free hydrogen atom sat- 
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FIG. 3: Quantum photo-absorption spectrum at the scaled 
energy E = —1.4 and the scaled electric field strength F = 
0.1. The initial state is |2p0), the light is polarized along 
the magnetic field axis. The plot shows the squared dipole 
matrix elements, which for graphical reasons are multiplied 
by it;. The strengths of the extraordinarily strong lines of the 
lowest n-manifolds at w < 7.5 are scaled down by a factor of 
0.2. 
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(57) 



so that in the scaled spectrum the unperturbed n- 
manifolds appear equidistantly spaced at 



2En 



(58) 



These spacings can clearly be recognized in figure |3j. At 
low values of w, neighboring n-manifolds are isolated. 
Furthermore, in this region the magnetic quantum num- 
ber to is nearly conserved. This is apparent from the fact 
that each n-manifold contains a central group of strong 
levels corresponding to to — 0, which can be excited even 
at F = 0, and adjacent groups of considerably weaker lev- 
els with to = ±1. Levels with higher magnetic quantum 
numbers are too weak in this region to be seen in the 
figure. At higher values of w, the conservation of m is 
violated, and individual n-manifolds acquire strong side 
bands. At even higher w, different n-manifolds strongly 
overlap. Throughout the spectral range shown, groups of 
strong lines indicating the centers of different n-manifolds 
are clearly discernible. 
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V. LOW-RESOLUTION SEMICLASSICAL 
SPECTRA 

A semiclassical approximation to a scaled photo- 
absorption spectru m is ob tained if the closed-orbit theory 
formulae of section IIIE| are rewritten in terms of scaled 
quantities, viz. 



g osc (w) — — Ac.o. cxp (iwS c 



with 



(59) 



(60) 



When low-resolution photo-absorption spectra are to 
be calculated from (^9|), a method of cut-off must be 
adopted to deal with the divergence of the semiclassical 
closed-orbit sum. For this section, we choose a Gaussian 
cut-off, i.e. (|59|) is replaced with 



9° sc 



(to) = — Ac.o. exp 



iwS c . . - 




(61) 



so that orbits with scaled actions larger than the cut- 
off action a are smoothly suppressed. This smoothing 
corresponds to a convolution of the quantum signal with 
a Gaussian of width X/a. 

To facilitate the comparison of (^lj) with the convo- 
luted quantum spectrum, we added the smooth part of 
the spectrum to <?° sc , which was calculated by convo- 
lving the quantum spectrum with a Gaussian of width 
1/ct = 1. This function is broad enough to wipe out the 
distinction between neighboring principal quantum num- 
bers. Results obtained for a = 20 and a = 50 are shown 
in figure [l] In both cases it is apparent that the large- 
scale structure of equidistant principal quantum numbers 
is well reproduced by the semiclassical approximation. In 
the quantum spectra, the substructure of the individual 
rt-shells can be discerned to a certain degree, given by 
the smoothing width 1/a. In the case of a = 20, much 
of this fine structure is also present in the semiclassical 
spectrum, but often the agreement is not good quanti- 
tatively. In particular, the peaks corresponding to the 
lowest n-manifolds are considerably wider in the semi- 
classical than in the quantum spectrum. 

If the cut-off action is increased to a — 50, finer details 
are resolved in the quantum spectrum. At the same time, 
the semiclassical closed-orbit sum becomes more oscilla- 
tory to reproduce this hne structure. It appears, how- 
ever, to be somewhat over-oscillatory, developing struc- 
tures absent from the quantum spectrum. This type of 
behavior is typical of closed-orbit sums in non-integrable 
systems. Thus, it can be questioned if the low-resolution 
closed-orbit sum can meaningfully be extended to even 
longer orbits. A high-resolution quantization based on 
the present semiclassical approximation will be presented 
in the following section. 
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FIG. 4: Smoothed quantum (upper halves) and semiclassical 
(lower halves, inverted) photo-absorption spectra with cut-off 
action (a) a = 20 and (b) a — 50. 



VI. HIGH-RESOLUTION SEMICLASSICAL 
SPECTRA 

For the calculation of a scaled semiclassical spectrum, 
the method of semiclassical qua ntization by harmonic in- 
version of 5 function signals |15|, ^6) can be applied. This 
technique requires the inclusion of closed orbits up to 
a maximum scaled action, i.e. it replaces the Gaussian 
cut-off used for the low-resolution semiclassical spectra 
presented in the previous section with a rectangular cut- 
off. A rough estimate for the required cut-off action can 
be obtained by means of perturbation theory fell . 



Stnax = —8-kEi 



i.e. w 



(62) 
15.06, this 



For the case E = —1.4 and n - 
estimate yields S' max /27r w 50. 

According to ([32]), to compute levels at high quan- 
tum numbers n a long semiclassical signal is needed, 
which can be hard or even impossible to obtain. We 
calculated closed orbits up to S max /2iT = 200, so that 
the orbital data is available for nearly 18,000 closed- 
orbit multiplcts. However, for reasons to be described 
in section VII a useful semiclassical signal can be con- 
structed up to S max /2ir w 60 — 70 only, so that, from 
the above estimate, the semiclassical calculation cannot 
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reach manifolds much higher than n = 10. On the other 
hand, the semiclassical approximation must be expected 
to yield more accurate results for higher quantum num- 
bers. Thus, when a high-resolution semiclassical spec- 
trum is to be calculated, a compromise must be made 
between the contradictory requirements of describing a 
spectral region at sufficiently high quantum numbers and 
with a sufficiently low spectral density. 

For the harmonic analysis of the closed-orbit sum we 
applied the method of 5 function decimated signal diag- 
onalization 26, p7|, which yields not only semiclassical 



eigenvalues and amplitudes, but also an error parameter 
estimating the precision of the eigenvalues. Results ob- 
tained for E = —1.4 and F — 0.1 with a signal length of 
<S m ax/27r = 60 are compiled in table The table contains 
the quantum eigenvalues of w and their dipole matrix ele- 
ments for levels satisfying (2pO\D\f)' 2 > 0.7. It is obvious 
at a glance that out of the multitude of spectral lines with 
intensities varying over many orders of magnitude (most 
of which are not contained in the table) only the strongest 
lines were detected in the semiclassical spectrum. The 
semiclassical eigenvalues given are characterized by hav- 
ing small imaginary parts, small error parameters and 
large amplitudes as well as being stable with respect to 
a variation of numerical parameters. The calculation op- 
erates at the edge of convergence, and in a few cases one 
can be in doubt whether a level should be included ac- 
cording to these fairly "soft" criteria, but in general a 
clear decision can be made. Semiclassical values for the 
transition strengths are not given because they are not 
reasonably well converged and depend strongly on the 
numerical parameters. 

One might expect that in each n-manifold it is the 
strongest lines that are detected semiclassically, and in 
general this expectation is confirmed by the numerical 
data. This can clearly be seen, e.g., in the manifold 
n = 6, which contains the most stably converged lines 
in the spectrum. There are, however, a few conspicu- 
ous exceptions, e.g. at n = 7, where strong lines are 
missing whereas comparatively weak lines are found. For 
n — 5, no lines at all can be computed from the given 
semiclassical signal. If the signal length is decreased to 
S , max /27r = 50, the three strongest lines appear in the 
spectrum in this manifold. 

At higher n, the number of strong lines in the quan- 
tum spectrum increases. So does the number of lines in 
the semiclassical spectrum until n = 11, where only three 
semiclassical lines are found. They appear rather arbi- 
trarily scattered across the quantum spectrum, and their 
convergence is notably worse than in lower manifolds. It 
is clear that in this n-shell the semiclassical quantization 
with the given signal is about to break down. At n = 12, 
no lines can be detected semiclassically. As, from the 
above discussion, this failure was to be expected because 
the required signal length becomes too large, the obvi- 
ous way to improve convergence seems to be to use a 
longer signal. However, if the signal length is increased 
to S max /2n = 70, no reasonably converged semiclassi- 
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FIG. 5: Stability determinants of vibrators as a function of 
the action for E = -1.4, F = 0.1. 



cal lines can be found in any n-manifold. Neither are 
results improved if the technique of harmonic inversion 
of cross-correlated closed-orbit sums [|l6[ Q is applied. 
This method has proven powerful in reducing the sig- 
nal length required in a semiclassical quantization. In 
the present case, however, because the cross-correlation 
increases the total number of frequencies obtained from 
the harmonic inversion, the true eigenvalues are hidden 
among a multitude of spurious frequencies, and no useful 
results can be obtained. 

For the time being, therefore, the results given in ta- 
ble H represent what can be achieved in the semiclassical 
quantization of the crossed-fields hydrogen atom. They 
confirm the applicability of the closed-orbit theory ap- 
proach in principle, but they also reveal a fundamental 
problem in its present formulation. From the analysis of 
the ideal test signal it is clear that the signal length avail- 
able is sufficient for a stable signal analysis. Thus, if the 
semiclassical results are not good, the semiclassical signal 
itself, rather than the signal analysis, must be to blame. 
This conclusion is confirmed by the observation that an 
increased signal length destroys the results rather than 
improves them. We therefore start searching for a flaw 
in the construction of the semiclassical photo-absorption 
spectrum. 

A conspicuous problem lies in the fact that the set 
of closed orbits available is incomplete. In no series of 
rotators or vibrators can arbitrarily long orbits be cal- 
culated. In the case of vanishing electric field there is 
a critical angle i? c which the starting angles of both ro- 
tators and vibrators approach as the orbits get longer. 
This convergence indicates that the orbits approach a 
separatrix between two families of tori in phase space. If 
sufficiently long orbits are studied, there are many closed 
orbits with very similar initial conditions, so that the 
numerical search for closed orbits must eventually fail. 

The region of phase space where the unknown orbits 
are located is lying close to a separatrix, so that it is 
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TABLE I: Semiclassical and quantum eigenvalues Wf of the scaling parameter for E — —1.4 and F - 
description. The dipole matrix elements {2p0jDj/} 2 were obtained from a quantum spectrum. 



0.1. See text for a detailed 



highly unstable. The orbits can therefore be expected 
not to contribute much to the semiclassical signal. The 
magnitude of an orbit's contribution to the closed-orbit 
sum ([33]) is determined mainly by its stability determi- 
nant M. Figure || shows the stability determinants of 
the vibrator orbits for E = —1.4, F = 0.1 as a func- 
tion of the scaled action. Different series of vibrators can 
clearly be discerned in the plot. It is indeed unstable or- 
bits with large M that are missing in the data set, but on 
the other hand the stability determinants of the missing 
orbits are not large enough to regard the corresponding 
semiclassical amplitudes as negligibly small. Because a 
vast majority of orbits has small M and was found, one 
can still hope that useful results can be obtained from 
the semiclassical signal, at least for quantum states not 
located in the separatrix region in phase space, but it is 



clear that the quality of the semiclassical signal is reduced 
by its incompleteness. 

To assess in detail the detrimental effect of the missing 
orbits and of any other sources of error that may ex- 
ist, we carry out a semiclassical analysis of the quantum 
spectrum. 



VII. SEMICLASSICAL RECURRENCE 
SPECTRA 

According to equation (p9|), in a scaled photo-absorp- 
tion spectrum every closed orbit contributes a purely si- 
nusoidal modulation to wg(w). This contribution can 
be extracted from the spectrum either by a conven- 
tional Fourier transform or by means of a high-resolution 
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method. The spectral analysis yields information about 
classical orbits returning to the nucleus. For this rea- 
son, the transformed spectrum is referred to as a recur- 
rence spectrum. High-resolution methods [H extract the 
scaled actions and scaled semiclassical amplitudes of in- 
dividual orbits and thus yield more complete information 
about the semiclassical spectrum than the Fourier trans- 
form, but they fail if the average density of closed orbits 
per unit of scaled action is too large. By contrast, due 
to its linearity the Fourier transform can be applied to 
any part of the recurrence spectrum with equal ease and 
numerical stability, irrespective of the spectral density. 
In dense regions, it will not be able to identify individual 
closed orbits, but it will nevertheless yield a recurrence 
spectrum that can be compared to the classical data. 
In this section we will present results obtained by both 
the Fourier transform and a high-resolution method. The 
semiclassical recurrence spectra will be compared to clas- 
sical results in order to identify the reason why the semi- 
classical signal is only partially suitable to a semiclassical 
quantization. 

Using either method, it is essential to note that the 
semiclassical closed-orbit formula cannot be expected to 
yield accurate results for the lowest levels. Thus, the 
low n-manifolds must be excluded from the semiclas- 
sical analysis, i.e. the analysis is based on the quan- 
tum spectrum given in an interval [w m i E ,w max ] instead 
of [0,w max ]. Furthermore, to minimize the impact of 
boundary effects due to the finite length of the semiclas- 
sical spectrum, a smooth Gaussian cut-off with width k 
centered at wq — (w m i n + w max )/2 is introduced. The 
smoothing replaces the peaks of the semiclassical recur- 
rence spectrum by Gaussians of width 1/k. The recur- 
rence spectra presented here were calculated from the 
quantum spectrum shown in figure for E = —1.4 and 
F = 0.1, with w min = 20, iu,nax = 100, and n — 10. For 
the high-resolution recurrence spectra, the method of S 
function decimated signal diagonalization was used. 

For low scaled actions, where only few closed orbits ex- 
ist, the high-resolution analysis can be applied. Results 
are shown in figure ^, which compares both the scaled 
actions and the real and imaginary parts of the semiclas- 
sical amplitudes extracted from the quantum spectrum 
to the classical results. For most closed orbits, the agree- 
ment is excellent. Exceptions occur for the shortest or- 
bits, where the actions of rotator and vibrator orbits are 
too similar to be resolved by the harmonic inversion. At 
somewhat larger actions, the three orbits in each group 
fall apart into two rotator orbits with similar actions and 
a vibrator orbit with a slightly larger action. 

These observations can be made even more clearly if 
the absolute values of the amplitudes are considered. 
They are shown in figure [?], where the results of the 
high-resolution analysis are also compared to those of the 
Fourier transform. Notice that for the Fourier transform 
the semiclassical amplitude is given by the area under a 
peak rather than the peak height, so that an immediate 
comparison to the high-resolution results is difficult. In 
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FIG. 6: High-resolution recurrence spectrum for E = —1.4 
and F = 0.1. Sticks and squares: semiclassical closed-orbit 
amplitudes, stars: harmonic inversion of the quantum spec- 
trum. 



figure [?], the Fourier transform was arbitrarily scaled so 
that the peak heights roughly match the values of the 
high-resolution amplitudes. For isolated orbits identified 
both in the Fourier transform and the high-resolution 
spectrum, the agreement between the two methods is ex- 
cellent. Where several peaks overlap in the semiclassical 
spectrum, no direct comparison is possible because the 
peak phases cannot be determined from the figure. 

Figure also extends the results shown in figure || to 
higher actions. In this region the density of closed or- 
bits starts to increase because, on the one hand, rotators 
of the first series exist and, on the other, bifurcations of 
closed orbits generate additional orbits. Apart from the 
fact that many orbits cannot be identified individually 
even by the high-resolution method, the most conspicu- 
ous feature of figure ^ is that for many orbits the semi- 
classical amplitudes calculated from the classical data are 
considerably larger than those extracted from the quan- 
tum spectrum. In some cases, this is obvious at a glance, 
but a closer inspection of the figure reveals that this phe- 
nomenon is rather common. Som e specific cases will be 
described in detail in section VIII. 



The occurrence of exceedingly large semiclassical am- 
plitudes is a well-known problem of both closed-orbit and 
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FIG. 7: Absolute value of the recurrence spectrum. Sticks 
and squares: semiclassical closed-orbit amplitudes, stars: har- 
monic inversion of the quantum spectrum. Solid curve, in- 
verted: Fourier transform (arbitrary units). 



periodic-orbit theory. It is associated with bifurcations 
of classical orbits, where, in the case of closed orbits, the 
stability determinant M vanishes and the closed-orbit 
amplitude (35) diverges. Close to the bifurcation, M 
is small. The semiclassical amplitude of the bifurcat- 
ing orbit is therefore large and exceeds the value deter- 
mined from the quantum spectrum. In a classical con- 
text, we have shown previously JT^ ] that vanishing M is 
both a necessary and sufficient condition for a bifurcation 
of closed orbits. In the context of semiclassical closed- 
orbit theory, it is necessary to overcome the divergence 
of the closed-orbit formula occurring close to a bifurca- 
tion. This problem will be addressed in section VIII, 



after the impact of the bifurcations on the semiclassical 
signal at hand has been investigated further. 

Whereas, in figure ^ the vibrator orbits are sufficiently 
isolated to be resolved by both the harmonic inversion 
and the Fourier transform across the entire range of ac- 
tions, the rotators occur in groups of several orbits having 
nearly identical actions. They are not resolved properly 
by either method. Instead, the Fourier transform pro- 
duces peaks describing the collective contribution of the 
orbits in a group. The harmonic inversion fits this contri- 
bution with fewer actions and amplitudes than the actual 





FIG. 8: Absolute value R(S) of the recurrence spectrum with 
k — 10 (see text). Upper part: Fourier transform of the quan- 
tum spectrum, lower part (inverted): smoothed semiclassical 
recurrence spectrum. 



number of orbits. Although the results can be expected 
to reproduce the quantum spectrum fairly well, the prin- 
cipal virtue of the high-resolution analysis - that it is 
capable of giving individual rather than collective contri- 
butions - is lost. It is therefore pointless to extend the 
high-resolution analysis to higher actions unless a signif- 
icantly longer quantum spectrum can be obtained, and 
only the Fourier transform will be used henceforth. 

Figure || displays the Fourier recurrence spectrum with 
smoothing n — 10 for scaled actions up to S/2ir = 100 
and compares it to the semiclassical spectrum. These re- 
sults extend the semiclassical analysis of quantum spec- 
tra to significantly longer orbits than investigated in pre- 
vious studies. They allow a verification of closed-orbit 
theory all the way up to the long orbits. It is immediately 
apparent from the figure that the quantum recurrence 
spectrum retains its pronounced peak structure. This is 
to be expected from closed-orbit theory, and indeed the 
peak locations are given by the actions of closed orbits 
for long as well as for short orbits. The basic idea of 
closed-orbit theory that recurrence peaks are related to 
classical closed orbits is therefore confirmed in principle 
even for very long orbits. 

Even for the largest actions considered, the quantum 
and semiclassical recurrence spectra agree quantitatively 
for some peaks. For most peaks, however, the peak 
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heights in the quantum and semiclassical spectra dis- 
agree. There are quantum peaks that are smaller in the 
semiclassical spectrum or even completely absent. They 
can be attributed to missing orbits. On the other hand, 
in many cases the semiclassical peaks are significantly 
higher than the quantum peaks, sometimes by several or- 
ders of magnitude. Exceedingly high peaks can be traced 
back to bifurcations of closed orbits if the possibility is ig- 
nored that a quantum peak can be small because orbits 
missing in the semiclassical spectrum interfere destruc- 
tively with the orbits present. This latter mechanism be- 
comes the more implausible the larger the semiclassical 
peak is in comparison to the quantum peak. 

Taken together, the effects of missing orbits and of bi- 
furcating orbits distort the semiclassical recurrence spec- 
trum to the point where it can no longer be expected 
to provide a suitable basis for a quantization. A close 
inspection of the recurrence spectrum suggests that the 
problem posed by bifurcating orbits is more severe. Ex- 
ceedingly high peaks do not only occur frequently, but in 
addition the very fact that they are high increases their 
detrimental effect on the semiclassical photo-absorption 
spectrum. Unless a suitable scheme for dealing with bi- 
furcating orbits can be devised, no improvement of the 
semiclassical signal can be expected. We therefore turn 
to a description of the semiclassical treatment of bifur- 
cating orbits by means of uniform approximations. 



VIII. UNIFORM APPROXIMATIONS 
A. The construction of uniform approximations 

Exceedingly large contributions of single orbits to a 
semiclassical spectrum arise when the orbits are too close 
to a bifurcation to be regarded as isolated, as is implicitly 
assumed by the stationary-phase approximation used in 
the derivation of the closed-orbit formula. Uniform ap- 
proximations furnish a collective contribution of all or- 
bits involved in a bifurcation. This solution was first 
suggested by Ozorio de Almeida and Hannay [2{| in the 
context of periodic-orbit theory. Their original approach 
was extended by different authors |3(], [H], |3^, [33| , so that 
today uniform approximations are a well-established tool 
of semiclassical physics. In reference jl7]], we identified 
two types of generic closed-orbit bifurcations of codimen- 
sion one. The pertinent uniform semiclassical approxima- 
tions will be derived in what follows. 

In most cases of interest, a bifurcation destroys real or- 
bits and turns them into complex ghost orbits that exist 
in the complexified classical phase space. Ghost orbits 
can yield palpable contributions to semiclassical spectra 
p3L M . In particular, their knowledge is essential for the 
construction of uniform approximation. For the generic 
closed-orbit bifurcations, the ghost orbits were described 
along with the real orbits in reference (T^| . 

Of particular importance is the observation that bi- 
furcations of codimension higher than one are relevant 



to semiclassics, although on a classical level they are 
not generically encountered. They appear as sequences 
of generic bifurcations, which, if the individual bifurca- 
tions are sufficiently close, must be described collectively 
by a single uniform approximation. Several examples 
of uniform approximations for these complicated bifur- 
cation scenarios have been described in the literature 

The principal requirement a uniform approximation 
must satisfy is to asymptotically reproduce the known 
isolated-orbits approximation when the distance from 
the bifurcation grows large, because in this limit the 
stationary-phase approximation can be expected to be 
accurate. In the following, we will advocate a some- 
what heuristic technique for the construction of a uni- 
form approximation, which is easy to handle and yields 
a smooth interpolation between the asymptotic isolated- 
orbits approximations on either side of the bifurcation. 
It will first be described in general terms. Subsequently, 
uniform approximations describing the generic types of 
codimension-one bifurcations of closed orbits will be de- 
rived. 

A bifurcation scenario is described by a normal form 
3> a (i) depending on n > 1 variables t and m > 1 param- 
eters a such that for any fixed value of the parameters a 
there are stationary points of & a (t) corresponding to the 
closed orbits involved in the bifurcation. The parameters 
a must then depend on the energy E to reproduce the 
bifurcations of the closed orbits. 

For the uniform approximation we make the ansatz 



with 



9(E) = 1(a) e tSo{E) 



1(a) = / d"tp(i)e l *» (t) . 



(63) 



(64) 



Here, the functions Sq(E) andp(i) as well as the parame- 
ter values a(E) have to be determined. All of them must 
be smooth functions of E. 

To find the asymptotic behavior of the uniform approx- 
imation ([33]) far from the bifurcations, (|64|) is evaluated 
in the stationary-phase approximation, which yields 



He) « £ 



(2 m y^ P (u) 

y/\Hess$ a (U)\ 



e »(S (E)+*o(tO) g-wrr/i/2 



(65) 

where the sum extends over all stationary points of 
$ a (t) that are real at the given a, Hess$ a is the Hes- 
sian determinant of $ a , and Vi is the number of negative 
eigenvalues of Hess $ Q (£i). This expression is supposed 
to reproduce the isolated-orbits approximation 



*(25)« Yl ME)' 



(66) 



In this case, the sum extends over all closed orbits in- 
volved in the bifurcation that are real at the given energy 
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E. If the normal form $ a (i) has been chosen suitably, 
there is a one-to-one correspondence between these or- 
bits and the stationary points fj. A comparison of J65| ) 
to ( |66| ) yields the conditions 



Si[E) = S (E) + 



and 



ME) 



{2mTl 2 p{U) 
v/lHcss*,,^)! 



—i-Kin/2 



(67) 



(68) 



These equations must be valid for real orbits. In most 
bifurcation scenarios, all orbits are real at least at certain 
energies. In these cases, it appears natural to postulate 
( |68| ) also to hold for ghost orbits. The parameter values 
one obtains are then smooth functions of the energy even 
at the bifurcations where the orbits become ghosts. In 
some instances, bifurcations involving only ghost orbits 
occur J36|, [38). In these cases, the condition (]6q ) still 
produces smoothly varying parameters and enforces the 
desired asymptotics. 

The numbers Vi of negative eigenvalues change dis- 
continuously at a bifurcation. For orbits which are real 
on either side of the bifurcation, so do the Maslov in- 
dices contained in the semiclassical amplitudes A%. These 
changes must compensate each other if the values p(ti) 
are to be continuous across the bifurcation. For these 
orbits, therefore, the change of Maslov index occurring 
in a bifurcation must be equal to the change in v± and 
can be determined from the normal form. For ghost or- 
bits, Maslov indices are not well defined classically. They 
must be chosen such as to make p(U) continuous. 

The normal form parameters a and the action Sq(E) 
can be determined from (|67|). They usually turn out to 
be unique. The amplitude function p(t), on the contrary, 
is unknown. Once the parameters a have been found, 
( |68| ) specifies its values p(ti) at the stationary points of 
$> a {t). These values, of course, do not suffice to identify 
p(t) uniquely, so that there is considerable freedom in the 
choice of p(t). Usually, if there are k orbits participating 
in the bifurcation scenario, we will approximate p(t) by 
a polynomial of degree fc — 1. This choice is justified by 
the observation that the uniform approximation is needed 
only close to a bifurcation, where all orbits are close to 
t = 0. Thus, in the spirit of the stationary-phase approx- 
imation, the dominant contributions to the integral (|64|) 
stem from the neighborhood of t = 0, whereas the regions 
of large t do not contribute. A suitable approximation 
to pit") must therefore be precise close to the origin. This 
is achieved by a Taylor series expansion, which leads to 
the polynomial ansatz. 

Simple as it might appear, however, this choice can 
bring about a mathematical difficulty: A polynomial p(t) 
diverges as t — ► oo, so that there is no guarantee that the 
integral flo4| ) will converge. If it does not, its divergence is 
an artefact of the choice of p(t), because by construction 
the regions of large t should not significantly influence 
the value of the integral. In this suitable regu- 

larization scheme must be applied. It can be justified 



by verifying that the regularized integral possesses the 
correct asymptotics. 

A slightly simpler form of the uniform approximation 
is obtained if the function p(t) is assumed to be a con- 
stant. This approximation does not exactly reproduce 
the desired asymptotics, but as the transition across the 
bifurcation mainly results in a change of the stationary 
points of Q a (t) rather than essential changes in p(t), it 
can be expected to capture the principal features. 

It is clear from the above description that there is a 
certain arbitrariness in the procedure. This arbitrariness 
can be reduced to the choice of a suitable amplitude func- 
tion p(t) , because by the splitting lemma and the classi- 
fication theorems of catastrophe theory |39| the uniform 
approximation can always be brought into the form ([63] ) 
by a suitable coordinate transformation, provided a nor- 
mal form & a (t) equivalent to the actual action function 
is given. 

In the following sections, uniform approximations for 
the two generic codimension-one bifurcations described 
in p7[ will be derived along the lines given here. They 
turn out to be analogous to those for isochronous and 
period-doubling bifurcations of periodic orbits given by 
Schomerus and Sieber [Mil. 



B. The fold catastrophe uniform approximation 

The simplest closed-orbit bifurcation is the creation of 
two orbits in a tangent bifurcation. It is described by the 
fold catastrophe 



$«(*) = ^ 3 - at . 



(69) 



This normal form has stationary points at t — ±\/a, 
which are real if a > 0. Its stationary values are (|70|) 



$(±V / ft) = T~ " 



3/2 



(70) 



By (p7|), the actions Si and S2 of the bifurcating orbits 
must satisfy 



51 = S (E)-^a 3 /\ 

5 2 = S Q (E) + ^a 3 / 2 . 



(71) 



For these equations to hold, one must assume Si < S2 
if the orbits are real and ImSi > 0, Im52 < if they 
are ghosts. These conditions determine how the orbits 
are to be associated with the stationary points of <& a (t) • 
Equation (fnl) can be solved for 



So(E) 



Si + S2 



and 



|o| = ( j \S 2 - Si 



2/3 



(72) 



(73) 
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The observation that the bifurcating orbits are real if 
a > and ghosts if a < fixes the sign of a. Both S {E) 
and a have thus be determined. 

For the semiclassical amplitudes, (^) yields 



mt/4 



With the ansatz 



(74) 



(75) 



for the amplitude function p(t) , we can solve for the pa- 
rameters po and pi to obtain 



Pa 



= V^\a\ 1/4 e-™/ 4 (A 1 +iA 2 ) , 



pi = Vtt - 



il/4 



-iir/A 



(.Ai - Ma) 



(76) 



The uniform approximation thus takes the form 

9(E) = ( Po Io+Pih)e lSo( - E) (77) 

with 



I k = — I dtt* e **-W . 
2tt 



(78) 



The integral Jo can be evaluated in terms of the Airy 
function as 

Jo - Ai(-a) , (79) 
whereas 7i is given by its derivative 

h = i-i-h = -«Ai'(-a) . (80) 
aa 

With these results, the uniform approximation ( |77| ) can 
be computed once the classical quantities S\,S2 and 
Ai,A 2 are known. After some rearrangements, ( f77| ) can 
be found to agree with the uniform approximation de- 
rived by Schomerus and Sieber |3l) for isochronous bi- 
furcations of periodic orbits, although its present form is 
much simpler. 



C. The cusp catastrophe uniform approximation 



action difference is given by the stationary value of <i> a (t) , 
which is a 2 /4, so that 



and 



— 'S'sym 'S'asym — ^G> 



a = ±2VAS . 



(82) 



(83) 



The parameter a has to be chosen positive if the asym- 
metric orbits are real, and negative otherwise. Here, AS" 
was assumed to be positive. If it is not, the normal form 
3> a (£) must be replaced with —<f> a (t), which changes the 
sign of the stationary values. 

Due to the reflection symmetry, the amplitude function 
must be an even function of t. We make the ansatz 



and solve 



P(t) =Po+P2t 2 
for the coefficients 

Po = 



(84) 



A p 1 

Z7T 



r/4 



Pi — ~Z i I -^asym V 2 iAsym 

The complete uniform approximation reads 
9(E) = J dtptye***® = poh+P2h 

with 

h = I dtt k e % *« [t) 



(85) 



(86) 



(87) 



The integral Iq can be evaluated analytically in terms of 
Bessel functions [Eol: 



T 71 /UT c -ia 2 /S Lttr/8 j a 

lo = 2 V N e [ e J -l/i I y 
+ sign a e 



J V4 -5- 



Although it is not apparent at first sight, I is a smooth 
function of a. This can be verified if the series expansion 



The normal form for the symmetrized cusp catastrophe 
is given by 



$«(«) = -i 4 - -at 2 
w 4 2 



(81) 



It has stationary points at t = and t = ±^/a and de- 
scribes a pitchfork bifurcation, where two asymmetric or- 
bits bifurcate off an orbit invariant under a reflection. We 
denote their actions and amplitudes by S sym , Sasym and 
-4 sym , „4 asym , respectively, where -4 a sym is understood to 
be the cumulative amplitude of both asymmetric orbits. 

As $ Q (< = 0) = 0, the reference action Sq(E) must be 
chosen equal to the action of the symmetric orbit. The 



J v (x) = (~Y r v (x) (89) 
with r v (x) a power series in x 2 is used. In terms of r v (x) , 



In = -e" 



**"r-v.(£) 



r 



t/8, 



■1/4 



(90) 



which is indeed smooth. The second integral I 2 can be 
evaluated from 
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da 



da 



-nr\/ a e 



2o 



= W8 J_i/ 4 (y)+signae— / 8 J 1/4 



J-5/4 I — I — J3/4 I — 



(91) 
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FIG. 9: Uniform approximation ( [si)) for the scaled spectrum 
at E = -1.4 and F = 0.2. 



This derivation contains an interchange of differentiation 
and integration which achieves a regularization of the 
divergent integral li. It can be justified by verifying that 
the asymptotic behavior of ( |9l| ) for a — > ±00 agrees with 
the stationary phase approximation to ( |S7| ) . 



IX. UNIFORMIZED RECURRENCE SPECTRA 

The formulae derived in the preceding sections give the 
uniform approximations directly in terms of the semiclas- 
sical actions and amplitudes. This circumstance makes 
them easy to apply to scaled spectra: we simply put 
S = wS and A c . . = w A c . .. As w is varied, the bi- 
furcation is not encountered because the classical me- 
chanics does not change, so that the isolated-orbits ap- 
proximation does not actually diverge. However, if w 
is small, the action differences between the bifurcating 
orbits are also small, so that the presence of the bifur- 
cation is felt and the isolated-orbits formula produces 
exceedingly large contributions. For large w, the action 
differences also grow large, so that the isolated-orbits ap- 
proximation should be recovered in the limit of large w. 

These findings are illustrated in figure O for a pitch- 
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FIG. 10: Contribution to the recurrence peak calculated from 
the uniform approximation (solid line) and the isolated-orbits 
approximation (dashed line) for the same bifurcation as in 
figure |, E = -1.4 and F = 0.1. 

fork bifurcation taking place in the first series of rotators 
at a repetition number /i — 57. At E = —1.4, the bi- 
furcation takes place at the scaled electric field strength 
F = 0.09014. The data shown in figure ^ was calculated 
for E — —1.4 and F = 0.2, which is sufficiently far away 
from the bifurcation for the asymptotic regieme to be 
reached within the range of w shown. As anticipated, in 
the limit of w — > 00 the complete uniform approximation 
agrees with the isolated-orbits formula. The simple ap- 
proximation also reproduces the beats correctly, but it 
has a smaller amplitude. 

The scaled uniform approximation can be used to im- 
prove the semiclassical recurrence spectrum, but this re- 
quires some effort: Whereas the isolated-orbits approx- 
imation yields 6 function peaks in the recurrence spec- 
trum, which are replaced with Gaussians due t o th e 
smoothing of the recurrence spectrum (see section |vTJ) , 
the uniform approximation is a complicated function of 
w. It must be subjected to a numerical Fourier trans- 
form in the same was as the quantum spectrum if its 
contribution to the recurrence spectrum is to be evalu- 
ated. Because a bifurcation involves orbits with roughly 
equal actions, the uniform approximation will produce a 
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recurrence peak at the appropriate action. An example is 
shown in figure |l(]. It was calculated for the bifurcation 
already describ ed i n figure 0. The Gaussian smoothing 
used in section VII was replaced with a rectangular win- 
dow, so that a number of side peaks appear. In this case, 
the Fourier transform of both the uniform approximation 
and the isolated-orbits approximation was taken over the 
rectangular window w £ [40,60]. The bifurcating orbits 
have the scaled action S/2w w 21.86, which is where the 
Fourier peaks are centered in both approximations. The 
peak produced by the uniform approximation is consid- 
erably smaller. 

If this uniformization procedure is carried out for all 
excessively high bifurcation peaks, it should be possible 
to bring the semiclassical recurrence spectrum in figure ^ 
into agreement with its quantum counterpart. In prac- 
tice, however, several obstacles stand in the way. First 
of all, in many cases ghost orbits must be included in 
the uniform approximation. They must be found and 
identified as pertinent to a given bifurcation before the 
uniformization can be performed. Furthermore, even if 
all relevant orbits are real, those orbits connected with 
each other in a bifurcation must be recognized in the data 
set. This is by no means an easy task. For example, if in 
a given series of rotators and for a given winding number 
a quartet of orbits appears, there are two different dou- 
blet orbits out of which they may have bifurcated, and it 
is not clear in general which of them must be taken for 
the uniform approximation. In a single case, this can be 
found out fairly comfortably by hand. If many orbits are 
to be classified, however, it is essential to do the group- 
ing automatically. We have not yet been able to devise 
a practical algorithm for this task, so that an automa- 
tized uniformization of all bifurcation peaks is presently 
impossible. 

Apart from these rather technical difficulties, there 
are also some obstacles of more fundamental impor- 
tance. Consider, e.g., the two high semiclassical peaks 
at S/2n w 25 in figure [?]. They are notably too high, and 
they are well-isolated from neighboring recurrence peaks, 
so that they may appear to be the ideal testing ground for 
the uniformization procedure. These peaks are generated 
by vibrators with repetition numbers fi = 41 and fi = 42, 
respectively. The pertinent bifurcation scenarios were 
described in figures 17 and 18 of (l7). The "simple" sce- 
nario taking place at fi — 41 consists of two orbits being 
generated in the rotational symmetry-breaking at F = 0, 
followed by a tangent bifurcation destroying one of them 
and a third orbit. To smooth this bifurcation peak, a 
uniform approximation describing the complete scenario 
must be found, which requires the construction of a suit- 
able normal form. Although a uniform approximation 
for the symmetry-breaking is available jll], 4^, 43 , the 
derivation of the pertinent normal form relies on prin- 
ciples different from the catastrophe theory classification 



used here, and it is not clear how these two can be united 
into a single normal form. Thus, the construction of a 
uniform approximation for this bifurcation scenario, and 
even more so for the more complicated scenario at /i = 42, 
remains an open problem to be solved in the future. It 
can be solved within the framework of uniformization 



presented in section VIII A , but will require a novel way 



of constructing normal forms. 

The approach to high-resolution semiclassical quanti- 
zation relies on the harmonic inversion of a Fourier trans- 
formed semiclassical spectrum, i.e. of a recurrence spec- 
trum. The above method of uniformizing the bifurcation- 
induced excessively high recurrence peaks in a semiclassi- 
cal spectrum would therefore, if it could be implemented 
systematically, also pave the way for the inclusion of uni- 
form approximations into a high-resolution semiclassical 
quantization, which has not been possible so far. We were 
able to demonstrate the feasibility of our method by way 
of example for the hydrogen atom in an electric field pl[ , 
which is less demanding classically. Its application to the 
crossed-ficlds hydrogen atom, however, remains open for 
future work. 



X. CONCLUSION 

For the first time, a high-resolution semiclassical quan- 
tization of the hydrogen atom in crossed electric and 
magnetic fields has been presented. It achieved the 
identification of the strong spectral lines in different n- 
manifolds. By means of a detailed semiclassical anal- 
ysis of the pertinent quantum spectrum, it was shown 
that bifurcations of closed orbits play a crucial role in 
the semiclassical spectrum and preclude the resolution 
of finer details in the semiclassical spectrum. They pose 
a particular challenge to the semiclassical quantization 
because they require a special treatment by uniform ap- 
proximations. 

A simple heuristic scheme for the construction of uni- 
form approximations has been proposed. Its simplicity 
and efficacy was demonstrated by a derivation of the 
uniform approximations for the codimension-one generic 
bifurcations of closed orbits. 

We have devised a general method for the inclusion of 
uniform approximations in a high-resolution semiclassical 
quantization by harmonic inversion. In a recent publica- 
tion [Q it was successfully applied to the hydrogen atom 
in an electric field. In the case of the crossed-fields hydro- 
gen atom, the diversity and complexity of the bifurcation 
scenarios encountered so far hinders the systematic im- 
plementation of the uniformization procedure. The treat- 
ment of all relevant bifurcations and the calculation of a 
detailed semiclassical spectrum thus remain challenging 
tasks for future studies. 
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